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Hybrid waves localized at hyperbolic metasurfaces 
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We reveal the existence of a new type of surface electromagnetic waves supported by hyperbolic 
metasurfaces, described by a conductivity tensor with an indefinite signature. We demonstrate that 
the spectrum of the hyperbolic metasurface waves consists of two branches corresponding to hybrid 
TE-TM waves with the polarization that varies from linear to elliptic or circular depending on the 
wave frequency and propagation direction. We analyze the effect of losses of the surface waves and 
derive the corresponding analytical asymptotic expressions. 


I. INTRODUCTION 

Metasurfaces are known as a two-dimensional ana¬ 
logue of metamaterials, and they offer unprecedented 
control over the light propagation, reflection, and refrac¬ 
tion EE- One of the main advantages of metasurfaces 
is that these structures are fully compatible with the 
modern planar fabrication technology and they can be 
readily integrated into the on-chip optical devices pre¬ 
serving the most of functionalities of three-dimensional 
metamaterials. In a general case, a metasurface can be 
described as a two-dimensional current characterized by 
a dispersive nonlocal two-dimensional conductivity ten¬ 
sor. At the same time, electromagnetic properties of a 
broad and constantly growing class of two-dimensional 
materials such as graphene, silicene, hexagonal boron ni¬ 
tride, can also be characterized by conductivity tensors. 
Thus, the physics of metasurfaces and optics of two- 
dimensional materials are tightly interconnected. Par¬ 
ticularly, it has been shown |3] that a graphene sheet can 
support transverse electric (TE) surface polaritons in the 
frequency region where the imaginary part of the surface 
conductivity becomes negative. Negative imaginary part 
of the conductivity corresponds to the capacitive surface 
impedance. At the same time, the existence of the TE 
surface waves at the capacitive impedance surfaces has 
been studied previously [4]. 

In this paper, we study a special class of metasurfaces 
characterized by a local diagonal anisotropic conductiv¬ 
ity tensor. Such metasurfaces can be regarded as a two- 
dimensional analogue of uniaxial crystals. Specifically, 
when the imaginary parts of the principal components 
of the conductivity tensors have different signs, a strong 
correspondence appears between these structures and hy¬ 
perbolic metamaterials [5]. Here we focus on the disper¬ 
sion and polarization properties of the localized waves 
supported by these metasurfaces. 
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It was shown that an anisotropic interface separating a 
hyperbolic metamaterial and vacuum can support a cer¬ 
tain class of plasmonic modes analogous to the Dyakonov 
surface states m- Dyakonov surface states [Id are lo¬ 
calized modes which can propagate in a narrow angle 
range along the interface of the anisotropic crystals. De¬ 
spite the theoretical prediction back in the 1980s, these 
modes have been experimentally demonstrated only re¬ 
cently m This is due to the fact that for the case of an 
interface of conventional anisotropic crystal these modes 
can propagate only in an extremely narrow range of an¬ 
gles, and thus it is hard to excite them experimentally. 
Nevertheless, these modes attract a significant scientific 
interest since they suggest a route for virtually lossless 
optical information transfer at the nanoscale, which is 
extremely important for the perspectives of on-chip op¬ 
tical data processing devices. Moreover, as was shown 
recently, the propagation direction of these modes can 
be effectively controlled by slight modification of the di¬ 
electric permittivities of the structure m- 

Here we show that, in a sharp contrast to bulk hy¬ 
perbolic metamaterials, the hyperbolic metasurfaces can 
support two types of surface modes at a single frequency. 
These modes originate from the coupling of the TM and 
TE polarized surface polaritons. A similar effect occurs 
in graphene sheets in the presence of a strong magnetic 
field perpendicular to the graphene layer [13] or on a 
metal substrate coated by a thin anisotropic dielectric 
film [14]. The corresponding surface waves have ellip¬ 
tic polarization which is essential for the construction of 
on-chip optical networks Earn 


II. MODEL 
A. Dispersion equation 

We consider two isotropic media with permittivities 
£i and 62 separated by an anisotropic non-chiral meta¬ 
surface [Fig. [lja)]. In the figure we suggest one of the 
possible design of hyperbolic metasurface representing a 
two-dimensional array of asymmetric subwavelength res¬ 
onators. Within the local homogenization approach the 
electromagnetic properties of such structure can be de¬ 
scribed by a two dimensional conductivity tensor: 
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FIG. 1. (a) Geometry of the problem. Upper [x > 0) half space with e\ and lower (x < 0) half space with £2 are separated 
by an anisotropic conducting layer. The principal axes of the conductivity tensor are shown by dashed lines. Surface wave 
propagates along z-direction. (b) Frequency dependence of imaginary parts of dimensionless conductivity tensor components 
<7j_ and cry. Parameters of conductivity tensor components are following: Q± — 1, Qy = 3, 7 = 0.05. 






(i) 


Here a± and cry are frequency dependent conductivities 
per unit length corresponding to the principal axes of the 
tensor. 

We will seek solution of the Maxwell’s equations in the 
form of a traveling wave propagating in the ^-direction 
and localized in the x-direction. The both electric and 
magnetic fields depend on the ^-coordinate and time t as 
exp (ik z z — iut). We assume that a is the angle between 
^-direction and one of principle axes of the tensor [see 
Fig. B a)]. Conductivity tensor is not diagonal in the 
chosen set of coordinates and can be written as 


and 0 and condition that electromagnetic field decays 
away from the interface 
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Here k\ 2 = k z — £ 1 , 2 <^ 2 /c 2 is inversed penetration depth 
of the surface wave into upper and lower medium. Sim¬ 
ilar equation describes dispersion of magnetoplasmons, 
surface waves in a two-dimensional electron gas in the 
presence of a strong DC magnetic field manansi. 


B. Conductivity tensor 


where 



( 2 ) 


(Jyy = ct_l cos 2 a + (j || sin 2 a ; (3) 

a zz = a l sin 2 a + cry cos 2 a ; (4) 

(Til — CF ±_ 

cr yz = cr zy = ---sm 2a. (5) 

Electric and magnetic fields (E and H) obey the fol¬ 
lowing boundary conditions on the metasurface: 

[n,H 2 ]-[n,H 1 ] = ^aE; (6) 

[n, E 2 ] — [n, EJ = 0. (7) 


Index 1 (2) corresponds to the upper (lower) half-space, 
n is a unit vector normal to the interface. 

Dispersion equation for surface waves can be obtained 
from Maxwell’s equations using boundary conditions 0 


Before analyzing dispersion equation (|8| let us discuss 
the form of conductivity tensor components a± and cry. 
Conductivity tensor describes a response of the layer on 
external electromagnetic field. In the common case, the 
response is not local in space and time that results in 
dependence of the conductivity on the frequency uj and 
wave vector k z . Equation (|8| was obtained using spa¬ 
tially local boundary conditions 00 - Therefore, direct 
substitution of a± and cry as function of k z into Eq. ([8| is 
not correct whereas Eq. 0 is still correct for frequency 
dependent a± and cry. 

Let us put forward the assumption that conductivity 
tensor components depend on the frequency according to 
the Drude-Lorentz model: 


ct s (uj) = A 
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Such dispersion shape is quite natural to many systems 
in optical, infrared, THz, and radio frequency ranges m ■ 


















3 


In Eq. © Q s are resonance frequencies along the prin¬ 
ciple axes of the conductivity tensor, 7 is a relaxation 
constant which is supposed to be isotropic. Constant A 
has a dimension of rad/s. Explicit expression of A is de¬ 
fined by the metasurface design. This constant can be 
excluded from the analysis with the help of the following 
dimensionless variables: 
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Real part of d_\_ and 07 is responsible for energy dissi¬ 
pation and imaginary part for the polarizability of the 
structure. Typical frequency dependence of Im(a^) and 
Im((T||) is shown in Fig.Jllb). One can see that signature 
of conductivity tensor ^ depends on the frequency. It is 
possible to distinguish three cases: (i) capacitive meta¬ 
surface when both Iul((t_l) and Im(5||) are negative; (ii) 
hyperbolic metasurface when Im(a ± )Im(5||) < 0; (iii) in¬ 
ductive metasurface when both Iui((t_l) and Im(d||) are 
positive. 
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FIG. 2. Dependence of k z on uj for the surface waves on 
hyperbolic metasurface for different propagation directions a. 
Two branches correspond to quasi-TM and quasi-TE surface 
plasmons. The inset shows the structure of dispersion curves 
at a « 90°. 


III. RESULTS AND DISCUSSION 


A. Dispersion of surface waves 


In order to analyze dispersion of surface waves which 
is described by Eq. © let us neglect dissipation of the 
energy in the system and put 7 = 0 . The case 7 7 ^ 0 is an¬ 


alyzed in Sec. IIIB For the sake of simplicity, further we 
will consider symmetric situation when £1 = £2 = £• As 
an example, let us consider the structure with resonance 
frequencies of conductivity tensor components U_l = 1 
and Cln =3. Dependence of wave vector k z on frequency 
uo and propagation direction a can be obtained analyt¬ 
ically from Eq. (| 8 |. This equation yields two solutions 
which correspond to hybrid polarized waves (quasi-TE 
and quasi-TM plasmons). Their dispersion for a = 60° 
is shown in Fig. [2] 

Surface waves of pure TE or TM polarization can prop¬ 
agate only along the principle axes of the conductivity 
tensor ( a = 0° and a = 90°). In this case, right part of 
Eq. ^ is equal to zero and the equation splits into two 
independent equations corresponding to two-dimensional 
TE and TM plasmons [21]. 

Frequency cutoff of quasi-TE plasmon is equal to zero 
and does not depend on a. Frequency cutoff uj c of quasi- 
TM plasmon belongs to the interval Q_\_ ^ uo c ^ fiy and 
depends on a. This dependence can be found from the 
equation: 


cot 2 
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Quasi TE-plasmon has a maximal frequency uj r at which 
it can propagate. Dependence of uo r on the a yields by 


the equation: 
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It follows from Eqs. (11) and ([ 12 ]) that simultaneous prop¬ 


agation of both surface waves at the same frequency is 
possible only if 
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This condition does not depend on the specific dispersion 
of <tj_ || and can be fulfiled for any hyperbolic metasurface 

m ■ ’ 

Dispersion curves of quasi-TM and TE plasmons for 
all angles are shown in Fig. [ 2 ] It follows from Eq. (J8| 
and can be seen from the figure that there is a frequency 
gap between to Uy at the angle a = 0°. The gap is 
squeezed as a tends to 7 t/2 . Structure of the dispersion 
curves at a close to 7 t/2 is shown in the inset of Fig. [ 2 J 
Presence of the anti-crossing means that there is a polar¬ 
ization switching along the dispersion curves. Therefore, 
notation ”quasi-TM” and ”quasi-TE” plasmons is just 
formal. 

Figure [ 2 ] contains full information about the dispersion 
of the surface waves on hyperbolic metasurface, but some 
peculiarities of the wave propagation associated, for ex¬ 
ample, with density of optical states or relative direction 
of phase and group velocities are best understood in k- 
space using isofrequency contours. Isofrequency contours 
for quasi-TE and quasi-TM plasmons are shown in Fig. 
[3] One can see that for quasi-TE plasmons, the contours 
have elliptical, 00 -shape or hyperbolic form depending 
on the frequency. Contours of equal frequency for quasi- 
TM plasmons have arc, rhombus, 8 -shaped or elliptical 































4 



FIG. 3. Contours of equal frequency on k-plane for quasi-TE (a) and quasi-TM (b) surface plasmons. Here k± and fey are 
component of the wave vector along the principle axes of the conductivity tensor. 


form depending on the frequency. The arc contours is 
observed for the hyperbolic regime when D_l < uj < fiy. 
End points of the arc correspond to the frequency cut¬ 
offs uj c which obey Eq. (11). Its solution represents 
fourth-order curve in k-plane, so-called hippopede [23] : 


uj 2 c = sin(<a ) 2 + fly cos(<a) 2 . (14) 

Discontinuity of the equal frequency contours at the finite 
points which takes place for the arcs in our case is unusual 
for bulk waves in three dimensional space but can be 
observed for the surface ones Ho]. Inside the hipopede 
quasi-TM plasmons are leaky-modes. 


B. Effect of losses 


Let us take into account energy dissipation and put 
7 / 0 in Eq. It results in finite propagation length 

of the surface waves which is proportional to 1 /Im (k z ). 
Frequency dependence of lm(k z ) for quasi-TE and TM 
plasmons is shown in Fig. [4j Using log-log scale makes 
it obvious that Im(^) ^ cu 4 at uj <C D_l and lm(k z ) ^ uj 
at uj D||. Dependence of Im {k z ) on uj can be obtained 
analytically from (| 8 | using perturbation theory. In the 
case of£< the losses can be written as 


Im(fc 2 )=^ 



(15) 


In the case of uj Dy the losses can be written as 


(16) 
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FIG. 4. Dependence of Im (k z ) on uj for quasi-TE and quasi- 
TM surface plasmons for a m 60° plotted on log-log scale. 
Dashed lines are asimptotics of the losses described. Curve A 
is described by Eq. (15); curve B is described by Eq. (16); 
curve C is described by Eq. (17). 


One can see that in contrast to the case of low fre¬ 
quencies, the losses do not depend on the propagation 
direction a and resonances frequencies and Dy. 

In the vicinity of the frequency cutoff uj c frequency 
dependence of lm(k z ) can be represented as 


lm(k z ) 


7 uj c o-(j( w c) sin 2 (a) + < 7 (w c ) cos 2 (a) 

2 1 +CT||(w c )5x(w c )/4 


■5u. (17) 


Here 5uj = uj — uj c and duo cj c . One can see from 


Im (k z ) = 27 a;. 
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FIG. 5. Polarization structure of hybrid surface waves on the 
hyperbolic metasurface. 


Eq. (17) that lm(k z ) 0 as Suj 0. It is a result 
of weak localization of quasi-TM plasmon near the fre¬ 
quency cutoff uj c . The opposite situation occurs for quasi- 
TE plasmon mode near uo r where this mode is strongly 
localized. 


C. Polarization properties 


Spatial inhomogeneity and partly longitudinal polar¬ 
ization of surface waves results in an unusual spatial 
distribution of their momentum and angular momentum 
density. In particular, surface waves can possess transver¬ 
sal Belinfante’s spin momentum [24} [25]. Anisotropy of 
the hyperbolic metasurface results in mixing of TE and 
TM surface waves and makes their polarization structure 
very manifold. 

It can be shown from the Maxwell’s equations that for 
the surface waves under consideration electric field com¬ 
ponents E z and E y are in-phase and E x has 7 r /2 phase 
delay. It means that electric field E rotates in the plane 
orthogonal to the metasurface so that the end of E draws 
an ellipse. In the common case, the plane of the ellipse 
is rotated on the angle (3 with respect to the wave vector 
k z (Fig. @). 

At low frequencies (5 D_l), there is only quasi-TE 

mode. Therefore, angle (3 is close to 90° and the ellipse 
is almost degenerated into a line segment. At high fre¬ 
quencies, where only quasi-TM mode propagates, angle 
f3 is near 0 ° and the ellipse represents a circle. 

Strong hybridization of quasi-TE and quasi-TM modes 
plasmon results in their unusual polarization. For exam¬ 
ple, for quasi-TM plasmon at a = 11.5° and uj = 2.95 
angle /3 = 87.3° and difference between semi-axis of the 
ellipse is less than 4%. Therefore, the wave has circular 
polarization. Absorption for this wave is quite low due 
to the vicinity of frequency cutoff uj c . In this case, fig¬ 
ure of merit (FOM) can be estimated using Eq. © or 
numerically 


FOM = Re( F } ss 1 • 10 5 for 7 = 0.05. 

Im(^) 


IV. CONCLUDING REMARKS 

We have presented a comprehensive analysis of surface 
waves propagating along hyperbolic metasurfaces. We 
have analysed dispersion, losses, and polarization prop¬ 
erties of such waves in the most general form, not speci¬ 
fying a specific design of the metasurface and describing 
its properties using the effective conductivity approach. 
Within this approach, the problem does not acquire a 
specific scale and, therefore, the results can be applied to 
different frequencies ranging from microwaves to ultravi¬ 
olet band. 

We have shown that the spectrum of waves supported 
by hyperbolic metasurfaces consists of two branches of 
hybrid TE-TM polarized modes, that can be classified as 
quasi-TE and quasi-TM plasmons. Dispersion properties 
of these waves are strongly anisotropic, and they have 
some similar features with magnetoplasmons and two- 
dimensional TE and TM plasmons. 

Analytical solutions of the problem allow detailed 
study of the surface wave properties. So, simple asymp¬ 
totic equations for the losses have been obtained near the 
frequency cutoff of quasi-TM plasmon mode and in the 
high and low frequency regions. Analysis of the equal 
frequency contours shows that their form and topology 
drastically depend on the frequency. The contours can 
have elliptic, hyperbolic, 8 -shaped, rhombus, or arc form. 
Multiplicity of equal frequency contours allows to fore¬ 
cast in hyperbolic metasurface such phenomena as neg¬ 
ative refraction mt self-collimation [27U29] . channeling 
of surface waves [501 EH], and large spontaneous emission 
enhancement of the quantum emitters due to the large 
density of states [321134] . 

We have shown that hyperbolic metasurfaces support 
simultaneous propagation of both quasi-TE and quasi- 
TM plasmon surface modes at the same frequency, and 
we derive the specific conditions for this to occur. Nei¬ 
ther in isotropic nor in chiral metasurfaces such a phe¬ 
nomenon is known Q3- Polarization structure of the sur¬ 
face waves can vary substantially so the polarization can 
change from linear to circular with different orientation 
of polarization plane. 

Unique electromagnetic properties of hyperbolic meta¬ 
surfaces make them quite promising for applications in 
many areas such as resonance sensing and detection, 
superlensing and near-field imaging, enhanced Raman 
spectroscopy, optical antennas, on-chip optical networks, 
etc. Taking into account their fabrication simplicity, rich 
functionality, and planar geometry it is possible to as¬ 
sert that hyperbolic metasurfaces can be a basis of many 
optical and optoelectronic devices. 
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